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ABSTRACT 

We perform several suites of highly detailed dynamical simulations to investigate 
the architectures of the 24 Sextantis and HD 200964 planetary systems. The 
best fit orbital solution for the two planets in the 24 Sex system places them on 
orbits with periods that lie very close to 2:1 commensurability, while that for 
the HD 200964 system places the two planets therein in orbits whose periods lie 
close to a 4:3 commensurability. In both cases, the proposed best-fit orbits are 
mutually crossing - a scenario that is only dynamically feasible if the planets are 
protected from close encounters by the effects of mutual mean motion resonance. 
Our simulations reveal that the best fit orbits for both systems lie within narrow 
islands of dynamical stability, and are surrounded by much larger regions of 
extreme instability. As such, we show that the planets are only feasible if they 
are currently trapped in mutual mean-motion resonance - the 2:1 resonance in 
the case of 24 Sex b and c, and the 4:3 resonance in the case of HD 200964 b and 
c. In both cases, the region of stability is strongest and most pronounced when 
the planetary orbits are mutually coplanar. As the inclination of planet c with 
respect to planet b is increased, the stability of both systems rapidly collapses. 

Subject headings: planetary systems - methods: numerical - planets and satel- 
lites: individual (24 Sex b, 24 Sex c, HD 200964b, HD 200964c - stars: individual 
(24 Sextantis, HD 200964 



1. Introduction 

Systems of multiple extrasolar planets provide rich laboratories for testing theories of 
planet formation and dynamical evolution. Radial- velocity planet-search programs, with ob- 
servational baselines approaching 20 years, often discover that apparent single-planet systems 
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show evidence for additional planets. These discoveries are due to a combination of increased 
time coverage and ongoing improvements in measurement precision arising from innovation 



and refinement in instru mentation, observing strategies, and analysis techniques ( jVogt et al. 



20101 ; iMayor et al.l 120111 ) . The former enables the detection of Jupiter-analog p 



orbital periods P > 10 yr (jWright et al.ll2009t IJones et al.ll2010l ; IWittenmver et al. 



anets with 



2012 ), and 



the latter permits the robust detection of planets of ever lower mass. I Wright et al.l (120091 ) 
showed that at least 28% of planetary systems surveyed by Doppler programs host multiple 
planets. Re cent results from the K epler spacecraft, with some 885 multiply-transiting planet 
candidates (IFabrycky et al.ll2012l ). have vastly increased the pool of these fascinating and 
valuable planetary systems. 

A number of recent studies have highlighted the need for observational detections of 
multiple-planet systems to be supported by dynamical investigations that test whether the 
orbits of the proposed planets are dynamically feasible. These investigations c an reveal that 
the p l anetary system configura t ion as reported is c atastrophically unstable ( 



2011 



Wittenmyer et al.l l2012al ; iHinse et al.l l2012al ; iGozdziewski et al.l 12012 



Horner et al. 



Horner et al. 



20121 ). Alternately, detailed dynamical simulations of systems close to mutual mean- motion 



resonance (M MR) can pro v ide important additiona l constr aints on the parameters of the 



planets (e.g. iHorner et al.l (120121 ); iRobertson et al.l (I2012al )). This is critical because any 



commensurability between the orbital periods of two planets can result in either extreme 
stability or instability, d epending on the precise orbits of the planets involved. For example, 
Robertson et al.l (2012a) pres ented revised orbits for the two Jovian planets in the HD 155358 
system ( Cochran et al. 2007h. potentially placin g them in mutual 2:1 MMR. The dynami- 
cal stability analysis of IRobertson et al.l ( l2012af ) showed that, for the orbital architecture 
presented in that work, the 2:1 MMR was indeed a region of long-term stability, serving as 
further evidence that the new orbital parameters are correct. In the HD 204313 three-planet 
system (jSegransan et al.ll2010l ). a similar dynamical analysis rev eals that the observed two 
outermost Jovian planets must be trapped in mutual 3:2 MMR (IRobertson et al.l 120121 ). an 
unusual architecture which dynamical mapping showed to be an island of extreme stability 
in the system. In that case, the dynamical analysis was critically necessary to constrain the 
HD 204313 system architecture, demonstrating that the two gas giant planets are locked in 
the 3:2 MMR. 

The intriguing results of these detailed stability investigations have prompted us to 



In recent years, a number of surveys (Hatzcs et al 



ems which ap 


pear 


to be in or near low 


(Hatzes et al. 


2005: 


Sato et al. 




2005: 



2006bl : iDollinger et al.ll2007l : I Johnson et al.ll2011bl : IWittenmver et al.ll2011bl ) have discovered 
a significant number of planetary systems orbiting intermediate-mass stars (M* > 1.5M Q ). 
These stars have proven to be a fertile hunting ground for interesting planetary systems. 
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Of particular interest are the two low-order resonan t giant planet pairs in the 24 Sex (2:1) 
and HD 200964 (4:3) systems (I Johnson et al.ll2011al ). For both systems, the best-fit radial- 
velocity solutions result in crossing orbits, an architecture which can only be dynamically 
stable on long timescales if the planets are protected from mutual close encounters by reso- 
nant motion. 



In their study of 24 Sex and HD 200964. iJohnson et al.l (j2011al ) performed small-scale re- 
body dynamical simulations which showed, as expected, that the long-term stable solutions 
were restricted to the 2:1 (24 Sex) and 4:3 (HD 200964) mean- motion resonances. The authors 
noted that their simulation results could not yet conclusively confirm that the planets were 
in resonance, and urged further, more detailed dynamical investigation. 

Resonant protection, such as that proposed for the 24 Sex and HD 200964 systems, is 
well known and studied in our Solar system. A wealth of Solar system objects move on orbits 



Hilda asteroids ( 


Franklin et al. 


1993 


), Jovian and Neptunian Trojans ( 


Morbidelli et al. 2005: 


Sheppard & Truiillo 


2006: 


Lvkawka & Horner 201G: 


lorner et al. 


20121) and. most famouslv. 


the dwarf planet Pluto and its brethren, the Plutinos 


Malhotra 


1995l:lFriedland 2001). Taken 



in concert with the growing catalogue of resonant exoplanets, these populations highlight 
the important role resonant dynamics plays in the formation and evolution of planetary 
systems. Indeed, astronomers studying the formation and evolution of our Solar system 
have learned a great deal about the extent, pace, and nature of the migration of the giant 



planet s through studies of the system's resona n t sma ll body population s (e.g. iMalhotra 
Jl995h : iMorbidelli et al.l (£3); iLvkawka et all fj2009h : broz et al.l (boilh ). revealing that 
the formation of our Solar system was most likely a relatively chaotic process. 

In this paper, we apply our highly detailed dynamical analysis techniques to the ex- 
tremely interesting 24 Sex and HD 200964 planetary systems. In Section 2, we briefly de- 
scribe the simulation methods and initial conditions. Section 3 gives the results, and we give 
our conclusions in Section 4. 



2. Numerical Methods 



To study the dynamics of the two planetary systems proposed in I Johnson et al.l (j2011al ). 
we performed two main suites of dynamical simulations - one for each of the planetary systems 
studied. We followed the strategy we have suc cessfully employed to s tudy the dynamics o f 



a number of other exopla n etary systems (e.g. iMarshall et al.l (120101 ) ; iHorner et al.l ( 1201 ll ) ; 



Wittenmyer et al.l (j2012al ): IHorner et al.l (120121 )). and followed the dynamical evolution of a 
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large number of different architectures f or each system us ing the Hybrid integrator within the 



n-body dynamics package MERCURY (jChamberslll999l ). In each case, we placed the better 



constrained of the two planets in question (24 Sex b and HD 200964 b) on its nominal best-fit 
orbit at the start of our integrations. The orbital parameters of the planets simulated in 
this work are shown in Table [T]0. For the other planet (24 Sex c and HD 200964 c), we then 
tested 41 unique orbital semi-major axes, distributed uniformly across the full ±3a range 
allowed by the uncertainties in that planet's orbit. For each of these 41 possible semi-major 
axes, we tested 41 unique eccentricities, again spread evenly across the full ±3cr range of 
allowed values. For each of these 1681 a — e values, we tested 15 values of the longitude of the 
planet's periastron, u, and 5 values of its mean anomaly, M, each spanning the appropriate 
±3<r error ranges. In this way, a total of 126,075 potential architectures were tested for each 
system. 

In each of our simulations, we followed the evolution of the two planets involved for a 
period of up to 100 Myr, until they were either ejected from their system, collided with one 
another, or were thrown into their central star. The times at which collisions and ejections 
occurred were recorded, which allowed us to create dynamical maps of the system's stability 
(Figures [T] and [3]) . 

In addition to the main suites of integrations discussed above, we performed five addi- 
tional suites of integrations for each of the two planetary systems to investigate the influence 
that the mutual inclinations of the planetary orbits would have on their stability. In this, we 



followed IWittenmyer et al.l (I2012al ). and considered cases where planet c was initially moving 
on an orbit inclined by 5, 15, 45, 135 and 180 degrees with respect to that of planet b. Due 
to the significant computational overhead in performing such runs, the resolution of these 
subsidiary investigations was lower than for the main runs. For each of these additional runs, 
a total of 11,025 trials were carried out, yielding the results shown in Figures 2 and 4. 



lr The orbital solutions for the planets were taken from the Exoplanets Data Explorer website, 
|http:/ /exoplanets.orgj on the 27th March 2012. We note that there are ve ry sma ll differences between 
the elements presented in that explorer and those presented in I Johnson et al.l (|2011a[) . These are, however, 
far smaller than the la errors on the values, and have no effect on our results or conclusions 
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3. Results and Discussion 

3.1. The 24 Sextantis System 

Figure [T] shows the results of our dynamical simulations for the 24 Sex 2:1 resonant 
system. At each point in the (a, e) grid, the small colored cell represents the mean survival 
time of 75 variants of the two-planet system - each with a unique initial combination of 
longitude of periastron and mean anomaly for the outer planet. The initial orbital parameters 
for the inner planet were held fixed, as noted in the previous section. The best-fit orbit for 
the outer planet is shown as an open box with crosshair s indicating the la un certainties in 



semi- major axis and eccentricity. The orbital solution of I Johnson et al.l (l2011aT ) lies directly 
on the narrow region of stable orbits, with mean survival times > 10 6 yr. This is strong 
evidence that the planets are truly in 2:1 resonance - that the resonance is required for the 
stability of the system. Nearly all of the surrounding parameter space is highly unstable, 
with survival times typically less than 10 4 yr. An additional small region of stability can be 
seen centered at a ~ 2.2 AU, for eccentricities below around 0.3. That region is the result 
of an overlapping web of weak, high-order resonances that congregate in t he region 2.174 - 



2.243 A U§. Interestingly, exactly the same feature can be seen in Figure 9 of iRobertson et al.l 



(j2012al ). The HD 155358 planetary system is somewhat analagous to that around 24 Sex, 
in that it features two planets that are most likely trapped in mutual 2:1 MMR. In that 
case, HD 155358c is most likely moving on an orbit that is somewhat less eccentric than that 
proposed for 24 Sex c, and so their Figure 9 shows the dynamical stability of that system 
to lower eccentricities. As the eccentricity of the outermost planet falls, the broad region 
of stability offered by these overlapped high-order resonances broadens until it merges with 
that offered by the protection of the 2:1 MMR. We note, however, that the stability region 
for 24 Sex c offered by those higher order resonances lies well away from the central ±lcr 
of the allowed orbital architectures for the system. Hence it seems far more reasonable to 
conclude that the planets in this system are, most likely, trapped in mutual 2:1 MMR. 

If the two planets were scattered to their p resent locations by a distant body, or by mu- 



tual c haotic interactions during their migration (IBarnes et all201ll : lMarzari fc Weidenschilling 



20021 ) . it is possible that they have some non-zero inclination relative to each other. To ex- 



2 More explicitly, the following commensurabilities are found in that region - 31:15 at 2.174 AU; 29:14 at 
2.177 AU; 27:13 at 2.181 AU; 25:12 at 2.186 AU; 23:11 at 2.191 AU; 21:10 at 2.197 AU; 19:9 at 2.205 AU; 17:8 
at 2.215 AU; 15:7 at 2.227 AU; 28:13 at 2.235 AU; 13:6 at 2.244 AU. At low eccentricities, such resonances 
can help to ensure the stability of orbits that would otherwi se be somewhat unstab l e, as has been obs e rved fo r 
both Solar system objects and exoplanetary systems (e.g. iLvkawka et al.l (|201 II ): IRobertson et al.l (|2012al) : 



Deck et al.l (120121) ). 
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plore the effect of mutually-inclined scenarios, we performed a subsidiary suite of integrations 
for a range of mutual inclination angles: 5, 15, 45, 135, and 180 degrees (i.e. coplanar but 
retrograde). These runs were set up as previously described, except at lower resolution: each 
scenario consists of a grid of 21 values of a, 21 of e, 5 of u, and 5 of M (11025 total trial 
systems). The results are given in Figure EJ and show that the system becomes generally 
more unstable when the planets depart from a prograde, coplanar configuration. Notably, 
the two retrograde scenarios had dramatically shorter lifetimes (panels e and f). We note 
that while the retrograde coplanar case (panel f ) shows a long-term stable region in the lower 
right, that region is quite far from the la uncertaint y on the orbit, much like our previous 
work on the proposed HU Aquarii planetary system (I Wittenmyer et al.ll2012al ). 



3.2. The HD 200964 System 



Figure [3] shows the results for the HD 200964 4:3 resonant system. As i n Figure [H the 



open b ox with crosshairs shows the best-fit parameters for the outer planet (IJohnson et al. 



2011al ). Once again, we see that the radial- velocity solution for this system, with the planets 
in a 4:3 resonance, lies within a narrow region of orbital stability surrounded by highly 
unstable parameter space. As was the case for the 24 Sex system, this dynamical map 
shows that the resonance is required for the system's stability. An interesting feature of the 
4:3 resonant protection is that the region of stability does not extend all the way to zero 
eccentricity - in other words, some small, non-zero eccentricity is required for HD 200964c 
to be dynamically stable within the 4:3 MMR with HD 200964b. This instability at very 
low eccentricities is observed in t he Solar system's Plutino population (trapped in 3:2 MMR 
with Neptune; see e.g. Fig. 6 of iRobertson et al.l (120121 )). and was also o bserved in some of 



the m ore extreme integrations of the planetary system orbiting HD 142 (IWittenmyer et al. 
2012h . 



At the far right-hand edge of the allowed range, for a > 2.05 AU, we see a highly stable 
region at all tested eccentricities. This feature is a common outcome of dynamical stability 
results at moderate and low eccentricities, loca ted just interior to the location of the 2:1 
MMR. A similar feature can be seen in Fig 1 of iHorner et al.l (120111 ). for the otherwise dy- 
namically unfeasible HU Aquarii planetary system. It typically represents the inner edge 
of the region for which dynamically stable solutions become the norm, rather than the ex- 
ception, apart from resonant interactions. At greater separations, this region of stability 
extends to ever greater eccentricities, since the boundary between stable and unstable so- 
lutions is determined, for the n on-resonant case, by the closest approach distance between 
the two planets in question. In IHorner et al.l (120111 ). this inner edge was discussed in terms 
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of the Hill Radius, Rh, of the more massive, innermost planet proposed in that work. The 
divide between stable and unstable orbits was found to follow a line that (roughly) followed 
a line of constant periastron distance for the outermost planet, centred on a periastron dis- 
tance between 3 and 5 Hill radii beyond the orbit of the innermost planet. For objects on 
near-circular orbits, the Hill Radius can be approximated as 

r h = a i/m 

where a and m are the semi-major axis and mass of the planet in question, and M is 
the mass of the central star. For HD 200964, the sharp boundary between unstable and 
stable orbits at around 2.05 AU once again lies between th ree and four Hill r adii b eyond the 



orbit of the innermost planet. Interestingly, we note that IChambers et al.l ( 119961 ) find that 
a system of two planets on low eccentricity, low inclination orbits "is stable with respect to 
close encounters if the initial semi-major axis difference, A, measured in mutual Hill radii, 
Rh, exceeds 2\/3, due to conservation of energy and angular momentum." The mutual Hill 
radius, Rhm is defined as 



p _ 3 / ( m b +m c \ ( a b +a c \ 

Rhm - y \~3M ){ 2 / 

where nib, rn c , o&, and m c are the masses and semi-major axes of planets b and c, and 
M is the mass of the central star. Given the criterion that the inner edge of the stable region 
should be found when the planets are separated by a distance of 2 a/3 times their mutual Hill 
radius, it is therefore trivial to work out where the inner edge of this stable region should 
be expected to lie, for the case where the orbits are circular. Holding the location of planet 
b fixed, we thus find that the inner edge of the stable region should lie at a = 2.15AU, a 
little more distant than that observed in our Figure [3j Once again, however, it should be 
noted that this stable region lies well beyond the central ±lcr region, and so represents a 
significantly less likely architecture for the HD 200964 system. 

As for the 24 Sex system, we also considered mutually-inclined scenarios, running a 
further series of simulations as described in Section 3.1. The results for the HD 200964 
system are the same: the mean lifetime decreases significantly when the planets are inclined 
by more than 90 degrees with respect to each other (retrograde orbits). Again, the most 
stable configuration was prograde and coplanar (panel a, identical to Figure 3). 



4. Conclusions 



Both the 24 Sex and HD 200964 systems host giant planets in close resonances (IJohnson et al 

2011aJ). We have performed detailed dynamical simulations testing the 3a range of allowed 
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parameter space for these two systems. Our results have further constrained the orbital 
parameters, with the best-fit solutions falling directly in narrow (~ la width) strips of long- 
term stability We also find that the stability of both systems is strongly dependent on the 
mutual inclinations of the planets involved, with coplanar orbits offering by far the greatest 
potential for dynamically stable solutions to be found. This work demonstrates the utility of 
such dynamical mapping for better understanding the architectures of multiple-planet sys- 
tems. The results of this work confirm that the resonant configurations are indeed required 
for long-term stability in the 24 Sex and HD 200964 systems. This adds them to a very short 
list of low-order-resonant exoplanetary systems, which are extremely valuable test cases for 
understanding giant-planet formation and migration processes. 
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Table 1. Planetary System Parameters 



Planet 


Period 
(days) 


To 

(JD-2400000) 


e 


(degrees) 


K 

(ms-i) 


M sin i 
(M Jup ) 


a 

(AU) 


24 Sex b 
24 Sex c 
HD 200964 b 
HD 200964 c 


455.2±3.2 
910±21 
613.8lJ;4 
825.0±3 * 


54758±30 
54941±30 
54916±30 
55029±130 


0.184±0.029 
0.412±0.064 

o.Mot8:8* 
o.i8il8;8S 


227±20 
172±9 
288l 4 [ 2 
182.6l^;I 


33.2±1.6 
23.5±2.9 
34.5l?i 
15.4±3.2 


1.6±0.2 
1.4±0.2 

1 - 84 -0.08 

0-895l8;^| 


1.338±0.024 
2.123±0.049 
1.601±0.002 
1.944±0.041 
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Fig. 1. — Dynamical stability for the 24 Sex system as a function of the initial semimajor 
axis and eccentricity of the outer planet. The nominal best-fit orbit for that planet is marked 
by the open square, and the la uncertainties are shown by the crosshairs. The 2:1 resonance 
appears as a narrow s trip of stability which coincides with the outer planet's best-fit orbit 
( 1 Johnson et al.ll2011ai ). 




Fig. 2. — Dynamical stability for the 24 Sex system, but for six values of inclination between 
the two planets. Panels (a) through (f) represent mutual inclinations of 0, 5, 15, 45, 135, 
and 180 degrees, respectively. Panel (a) is a duplicate of Figure 1, shown here for ease of 
comparison. As in previous figures, the color bar represents the log of the mean survival 
time. 
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Fig. 3. — Dynamical stability for the HD 200964 system as a function of the initial semimajor 
axis and eccentricity of the outer planet. The nominal best-fit orbit for that planet is marked 
by the open square, and the la uncertainties are shown by the crosshairs. The 4:3 resonance 
appears as a narrow s trip of stability which coincides with the outer planet's best-fit orbit 
( 1 Johnson et al.ll201fai ). 
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Fig. 4. — Same as Figure [2J but for the HD 200964 system. As for 24 Sex, mutual- inclination 
scenarios are much less stable than the prograde-coplanar scenario in panel (a). 



